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Abstract 

We have presented a model in which the new generalized Chaplygin gas interacts with 
matter. We find that there exists a stable scaling solution at late times in the evolution of 
the universe. Moreover, the phantom crossing scenario is observed in this model. 
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1 Introduction 



Astrophysical observations show that more then seventy percent of the cosmic energy 
density is contained in an unknown 'dark' sector commonly termed as 'dark energy' pQ 
(see [2] for reviews on dark energy). The remaining part of the total energy density is 
due to matter, which is also mostly dark [3J. This mysterious dark energy is specified by 
an equation of state (EoS) pd e = uJxPde, where pd e and pd e are the pressure and energy 
density of dark energy, while ujx is the corresponding dimensionless EoS parameter. In 
the presence of dark energy, the fabric of spacetime expands in an accelerated manner, 
implying pd e < and consequently u>x < 0. Since general relativity cannot satisfactorily 
explain the cosmic accelerated expansion, it motivates theorists either modifying curvature 
or the matter part in the Einstein field equations. In recent years, such options are 
carefully investigated in literature. Alternative gravity theories based on the modification 
of Einstein-Hilbert action include f(R) gravity [3], scalar-tensor gravity [5] and Lovelock 
gravity [Bj, to name a few. Similarly there are models in which only the matter term is 
modified, examples are Cardassian model [7j, the bulk viscous stress [8] and the anisotropic 
stress |9J. However all such modified models have several drawbacks in explaining the 
observational data. There is also a possibility that the dark energy admits an exotic EoS 
that manifests the observed accelerated expansion. We consider such a possibility by 
introducing an EoS based on the Chaplygin gas (CG). In the context of cosmology, the 
Chaplygin gas was first introduced by Kamenshchik et al [TO] . It is specified by p = —Lj p, 
where L is a constant. Its density evolution is given by 



where M is constant of integration and z is the redshift parameter. Interest in CG arose 
when it appeared that it gives a unified picture of dark energy and dark matter i.e. under 
certain constraints on L and M, expression (1) gives density evolution of matter at high 
redshifts and dark energy at low redshifts [11] . Other successes of CG is that it explains 
the recent phantom divide crossing [12], is consistent with the data of type la supernova 
[13] and the cosmic microwave background [14]. The CG emerges as an effective fluid 
associated with <i-branes [15] and can also be obtained from the Born-Infeld action [16] . 
Since matter and dark energy are the dominant components of the cosmic composition, 
it is natural to expect their mutual interaction at some scale. The exact nature of this 
interaction is still unexplained and the interaction may not necessarily be gravitational 
either. Cosmological models based on the interaction between dark energy and matter 
are termed 'interacting dark energy' in literature and are under thorough investigation 
[17j . In these models, either cosmic specie decays into the other depending on the sign 
of the coupling parameter involved. Recent interest in interacting dark energy models is 
also triggered from the astrophysical observations which show that the energy densities of 
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matter and dark energy are of the same order of magnitude i.e. r m = p m J pd e — 1. It leads 
to the 'cosmic coincidence problem' which asks the explanation of r m ~ 1 at the present 
time. Alternatively, why the dark energy parameter ujx is close to —1 in recent times. 
In the interacting dark energy scenario, there has been successful attempts in resolving 
this problem and stable attractor solutions of the Friedmann-Robertson- Walker (FRW) 
equations are obtained which give u>x — —1 closer to present time [18J. It is shown in 
P3] that the cubic corrections to the Hubble law, measured by distant supernovae type 
la, probes this interaction. Moreover, this interaction is controlled by third and higher 
derivatives of the scale factor. Moreover, cosmic microwave background observations lead 
to a constraint on the coupling parameter of the interaction as 0.05 < c < 0.2 [20] . 
Observationally the Abell cluster A586 provides evidence of the interaction between dark 
matter and dark energy |21j . 

In the context of field theory and particle physics, it is customary and appealing to 
interpret the dark energy as some sort of particles that interact with the particles of the 
standard model very weakly. The weakness of the interaction is required since dark energy 
particles have not been produced in the accelerators and because dark energy has not yet 
been decayed into lighter or massless fields such as photons. The interaction between dark 
energy and other particles cannot be arbitrary since this interaction gives a fifth force with 
a range A ~ 1/rriQ, where mo is the the mass of dark energy particle. It has been shown in 
[22] that an equation of state with u>d e < — 1 can be a signal that dark energy will decay in 
the future and the universe will stop accelerating. This conclusion is based in interpreting 
a u>de < — 1 as a signal of dark energy interaction with another fluid. In another paper [23] . 
it is shown that the mass of dark energy particle could be of the order ~ (10~ 33 eV) 2 . 
In the same study, it is proposed that the phantom particle can decay into one or more 
phantom plus an ordinary baryonic particle. Moreover, an ordinary particle may decay 
into phantoms plus other ordinary particles with a larger effective mass than the original. 
Thus the above discussion shows that model of interacting dark energy is supported by 
both theoretical arguments and observational evidences. 

Zhang and Zhu [23] used the Chaplygin gas in the interacting dark energy model and 
obtained the stable scaling solution of the FRW equations. Later on, their work was 
extended by Wu and Yu [25] for the generalized Chaplygin gas. We here extend these 
earlier studies by using the new generalized Chaplygin gas. 
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2 Modeling of dynamical system 



We start by assuming the background to be spatially homogeneous and isotropic FRW 
spacetime 

ds 2 = -dt 2 + a 2 {t)[dr 2 + r 2 {d0 2 + sin 2 6d<f>% (2) 

The equations of motion corresponding to FRW spacetime filled with the two component 
fluid are 

K 2 

H = —(Pncg + Pncg + Pm), (3) 

K 2 

H 2 = —{Pncg + Pm). (4) 

Here k 2 = 8tcG is the Einstein's gravitational constant and H = H(t) = a/a is the 
Hubble parameter. In this paper, we solve the FRW equations using the 'new generalized 
Chaplygin gas' (NCG) state equation, proposed by Zhang et al [26J. It is an extended 
form of the generalized Chaplygin gas and hence dubbed with the 'new'. The NCG 
model is dual to an interacting XCDM parameterization scenario, in which the interaction 
is determined by the parameter a. Here the X part corresponds to the quintessence 
{tjj x < —1/3), following the notation used in [26]. Since the observational data favors u>x 
to be in the range ( — 1.46, —0.78) [1, 2J, it has motivated to generalize the Chaplygin gas 
EoS to the NCG form to incorporate any X-type dark energy in the universe. 

p nC g = -^, A{a) = -u x Aa-^^\ (5) 

Pncg 

The density evolution of NCG is given by 

Pncg = [Aa-^ 1+ ^ {1+ ^ + Ba- 3 ( 1+a )]^, (6) 

where B is the constant of integration. The energy conservation equation for the dynam- 
ical system under consideration is 

Pncg + Pm + 3# (p m + p ncg + p ncg ) = 0. (7) 

Due to interaction between the two components, the energy conservation would not hold 
for the individual components, therefore the above conservation equation will break into 
two non- conserving equations: 

p ncg + 3H(p ncg + p ncg ) = -Q, (8) 
p m + 3Hp m = Q. (9) 

Here Q is the energy exchange term which is to be specified ad hoc. However from 
the dimensional considerations, it is obvious that Q should have dimensions of density 
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into inverse time, the later being chosen to be Hubble parameter. Thus we expect that 
Q = Q(Hp ncg , Hp m ) which upon expanding about densities in a Taylor series yields 
Q ~ H(p ncg + p m ) [27 \. We also insert a coupling parameter c in Q to determine the 
strength of the interaction, thus we have 

Q = 3Hc(p ncg + p m ). (10) 

From the observational data of 182 Gold type la supernova samples, CMB data from the 
three year WMAP survey and the baryonic acoustic oscillations from the Sloan Digital 
Sky Survey, it is estimated that the coupling parameter between dark matter and dark 
energy must be a small positive value (of the order unity), which satisfies the requirement 
for solving the cosmic coincidence problem and the second law of thermodynamics [28] . 
The positive c implies that the energy will flow from the NCG into matter. To study the 
dynamics of our system, we proceed by setting 

x = In a = -ln(l + z), (11) 

which is termed as the e-folding time parameter and z is the redshift parameter. Moreover, 
the density and pressure of NCG can be expressed by dimensionless variables u and v as 

= IW = ^ = ^, „ A. (12) 



The EoS parameter u>x is conventionally defined by 

Pncg 
Pncg 



uo x {x) = (13) 



which after using (12), becomes 

lu x (x) = -. (14) 
u 

The density parameters of NCG and dark matter are related as 

= = 1 - Qncg = 1 ~ U. (15) 

Using Eqs. (12), (13) and (15) in (8) and (9), we obtain 
du 

— = — 3c — 3v + 3uv, (16) 
ax 

dv Oil) , , . . . , , . , 

= 3— (u + v + c ) + 3v (1 + at )(1 + uj x ) + 3v (1 + v . 17 
ax u 

The critical points of the above system are obtained by equating Eqs. (16) and (17) to 
zero. The only critical point of the system is 

2-c + ux ( . 

2 + ujx 

v c = -(2 + lux). (19) 
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Notice that for ojx — — 1, our results reduce to those of [25] for the interacting generalized 
Chaplygin gas. Since in a spatially flat universe, the meaningful range is < u < 1, 
consequently < u c < 1 implying < c < 2 + ujx- Since uix < 0, therefore c > 1 can 
arise in certain cases. Note that the acceleration in the late evolution of the universe 
arises when 



aH 2 



< -1, 



Using Eqs. (12) to (15) in (20), one can write 



-(l + 3w). 



Using (19) in (22), we get 



Q 



3(2 + u x )). 



(20) 



(21) 



(22) 



For —1.6 < ujx < 0, we notice that q < 0, thus giving accelerated expansion of the 
universe (see Fig. 6). We further check the stability of the dynamical system (Eqs. 16 
and 17) about the critical point (u c ,v c ). To do this, we linearize the governing equations 
about the critical point i.e. u = u c + 5u and v = v c + 8v, we obtain 

d8u 



dx 
d5v 

dx 



3v c Su + 3(— 1 + u c )8v, 
—3a 



-(v c +cv c ) 



5u 



+ 



6av c 3ac 

3(1 + 2v c + a) + + + 3(1 + a){\ + u x ) 



U c U c 

The eigenvalues corresponding to the linearized system are 

1 



5v. 



(23) 



(24) 



A 



2(-2 + c 
-3. 



■[-6c(2 + uj x ) + 3(2 + a) (2 + u x ) 



'a{2 
1 



2(-2 + c-lu x 



w x ) 2 (-4c 2 + 4c(2 + uj x ) + a(2 + oo x ) 
-[-6c(2 + ux) + 3(2 + a) (2 + cu x ) 



+3^a(2 + 



-Ac 2 + 4c(2 + lux) + a(2 + u x ) 



(25) 



(26) 



To obtain stable critical point, the real parts of the eigenvalues must be negative. We 
notice that fora>0,0<c<2 + ojx and uj x < 0, the two eigenvalues A and \i will be 
negative, giving a stable attractor solution at (u c , v c ) in the late time evolution governed 
by the FRW equations. For instance, taking c = 0.3, ojx = —1.2 and a = 0.2 (generalized 
phantom energy), we have A = —1.86 and \i = —3.69; whilst for a = 1 (phantom energy), 
we have A = —1.64 and fi = —6.99. The negativity of the eigenvalues can also be proven 
as 



- 6c(2 + oj x ) + 3(2 + a)(2 + u x f = 3(2 + u x )(-2c + (2 + a) (2 + u x )) > 0, (27) 
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since a>0,0<c<2 + uox and ujx < 0. Now notice that both the eigenvalues will be 
negative if 



(-6c(2 + uj x ) + 3(2 + a) (2 + u x f > 9a(2 + cu x ) 2 (-4c 2 + 4c(2 + + a(2 + w x ) 2 ).(28) 



Therefore both the eigenvalues are indeed negative. In Figures 1 and 2, we have provided 
a pictorial relationship between u and t> for suitable choices of model parameters and 
initial conditions. It is interesting to notice that solutions converge to the same state 
under the given conditions. Figures 3 and 4 show that the functions u and v are stable 
solutions of the system i.e. the corresponding curves become flat at approximately x ~ 3 
and become globally flat afterwards. Initially these functions are unstable for small values 
of e-folding time x but in the late times, these become globally stable. It can also be seen 
that NCG state parameter ojx = v/u crosses ujx = — 1 at present time (x = 0) and enters 
the phantom regime u>x < —1 for large e-folding time (see Fig.5). 

3 Conclusion 

Model of interacting dark energy possesses enormous potential in resolving or at least sim- 
plifying several cosmological problems including the coincidence problem and the phantom 
crossing. In this model, these problems are analyzed by considering two major cosmic 
components namely dark matter and dark energy and assuming an interaction between 
them. Two non-conserving equations are derived and then non-dimensionalized. Due to 
non-linearity in the governing equations, we sought for numerical solution of the system. 
The stability of the solution is determined by performing stability analysis. In this paper, 
we have studied an interacting dark energy model dealing with the interaction between 
new generalized Chaplygin gas and the matter. The analysis is performed by doing the 
stability check of the dynamical system. We found that the dynamical system is stable 
about the only critical point of the system. This shows that stable stationary attractor 
solution exists in the late time evolution of the universe when the later one enters a steady 
state. Finally, this paper presents an extension of earlier work in |25j . 
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Figure 1: The phase diagram of the interacting new generalized Chaplygin gas model with 
c = 0.7, a = 0.07 and lox = —1-3. The curves correspond to the initial conditions u(— 2) = 
1.2,u(-2) = -0.2 (green); u(-2) = 1.3,«(-2) = -0.3 (blue); u(-2) = lA,v(-2) = -0.4 (red); 
u{-2) = 1.5 ! u(-2) = -0.5 (brown). 
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ure 2: The phase diagram of the interacting new generalized Chaplygin gas model with 
0.003, a = 0.3, ux = —1/3. The curves correspond to the initial conditions as given in Fig. 
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Figure 3: The parametric functions u and v are plotted against the e-folding time parameter 
x. The parameters are fixed as c = 0.3, a = 0.4 and ujx = — 1. The initial condition is 
u(0) = 2,v(0) = -2. 
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Figure 4: The parametric functions u and v are plotted against the e-folding time parameter 
x. The parameters are fixed as c = 0.5, a = 0.08 and ujx = —1/3. The initial condition is 
u(0) = 3,v(0) = -1. 
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Figure 5: The EoS parameter lux = v/u is plotted against x. The parameters are fixed as 
c = 0.5, a = 0.08 and lux = — 1- The initial conditions are n(0) = 1, u(0) = —1 (solid line); 
u(0) = 1.1, v(0) = -1.1 (dots); u(0) = 1.2, v(0) = -1.2 (dashes); n(0) = 1.3, v(0) = -1.3 (dash 
dot). 
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Figure 6: The deceleration parameter q is plotted against x. 
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